One-dimensional steady-state structures at relativistic interaction of laser radiation 
with overdense plasma for finite electron temperature 



OO 
O 

o 

» 

o 
O 

m 

Oh- 



, 

• ^ ■ 

^ : 

Oh. 



> 

in 

(N 



OO 

O 



X 



A. V. Korzhimanov^ and A. V. Kim^ 

'institute of Applied Physics RAS, 603950 Nizhny Novgorod, Russia 

One-dimensional steady-state plasma-field structures in overdense plasma are studied assuming 
that the electron temperature is uniform over plasma bulk and the ions are stationary. It is shown 
that there may exist solutions for electron distributions with cavitation regions in plasma under the 
action of ponderomotive force 
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The interaction of superstrong laser radiation with 
plasma has been recently arousing keen interest for sev- 
eral reasons Firstly, with the advance of laser tech- 
nologies radiation intensities up to 10^^ — 10^^ W/cm^ 
have been attained 0]- Secondly, wide opportunities 
open up for using new technologies for practical purposes, 
such as adronotherapy based on laser accelerated ions [1] , 
inertial nuclear fusion [F, ^ , generation of coherent X-ray 
and gamma-spectrum radiation , attosecond pulse gen- 
eration 0], and others. Finally, this field of research is 
quite complicated and poorly studied yet, hence there are 
still many problems to be solved, including fundamental 
ones. 

It is well known that to provide optical radiation at 
intensities of about 10^^ W/cm^ the electron motion 
in the wave field becomes essentially relativistic, which 
leads to pronounced changes in the dynamics of interac- 
tion of such radiation with plasma. In particular, owing 
to Lorentz mass increase the effective plasma frequency 
LUp = {ine^Ne/rriey/^ (Ne,me are the electron concen- 
tration and mass, respectively) grows and the plasma 
that was initially nontransparent at this frequency may 
become transparent iQj. Besides, the role of the magnetic 
component of the Lorentz force acting on the electrons 
from the direction of the electromagnetic wave increases 
considerably in the relativistic regime because this force 
is proportional to the Ve/c ratio, where Ve is electron 
velocity. This force evokes ponderomotive action of ra- 
diation on plasma which, at so high intensities, plays a 
decisive role in plasma dynamics upon the whole. 

As was mentioned above, the interaction of super-high 
laser radiation with plasma is a very complicated prob- 
lem and the main instrument for theoretical research 
in this field is numerical simulation on super-powerful 
computers. However, interpretation of results of these 
studies and problem formulation demand physical intu- 
ition based on simple analytical calculations. Therefore, 
of particular significance in terms of understanding the 
physics of the processes are the results obtained analyti- 
cally in some model problems. 

In the current work we study one-dimensional steady- 
state structures arising in overdense plasma normally 
irradiated by an intense relativistic laser pulse. This 
research is a logical continuation of the work done in 
[13, [m [H, [H, [iJI , in which analogous studies were car- 



ried out in the "cold" plasma approximation when the 
energy of electron thermal motion was negligible com- 
pared to the energy of their oscillatory motion in wave 
field. Within the above approximation the plasma-field 
structures were discontinuous as the small parameter of 
the higher derivative was neglected. In the present work 
we introduce this small parameter characterizing electron 
temperature, which enables constructing structures anal- 
ogous to the obtained earlier but sufficiently smooth. 

In our analysis we will use a joint and self-consistent 
system of Maxwell's equations for e.m. field and relativis- 
tic hydrodynamic equations for electron plasma compo- 
nent. Electron motion will be neglected under the as- 
sumption that the ion mass is much higher than the elec- 
tron mass, so that the times of their reactions are negli- 
gible compared to the intervals of the time of interaction 
of interest to us. For description of e.m. field we will 
introduce vector A and scalar ip potentials and will use 
the Coulomb gauge divA = 0. Then, in the approxima- 
tion of a steady-state monochromatic circularly polarized 
plane wave the vector potential may be represented in the 
following dimensionless form 



eA(z, t) 
moc^ 



(1) 



where Xo,yo are unit vectors of the x- and y-axes, z is 
the axis along which the wave is propagating, lu is wave 
frequency, and mo is electron rest mass. As was shown 
in |14l |. Maxwell's equations in plasma in this case may 
be reduced to the following equations 



dC^ 



a = 



= no{n - ZiUi), 



(2) 



(3) 



where C, = ujz/c; (f> = et^/moc^; s is the normalized e.m. 
energy flux density that is the parameter of the problem; 
no = Ane'^No/moLLj'^ is the plasma overdense parameter 
equal to the ratio of a certain characteristic initial elec- 
tron concentration A^o to its critical value at a given radi- 
ation frequency in the linear approximation (hereinafter 
we will consider only homogeneous plasma layers in which 
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Nq is taken to be the value of unperturbed electron con- 
centration); Zi is the ion charge (all the ions are sup- 
posed to have the same charge); n = Ns/Nq, rii = Ni/No 
are the normalized electron and ion concentrations, re- 
spectively; and 7 — VT+(? is electron relativistic factor 
(i.e., the ratio of their energy to rest energy; thermal en- 
ergy is neglected here). Thus, the above simplification 
enabled us to reduce Maxwell's equations to a nonlinear 
Helmholtz equation for vector potential complex ampli- 
tude modulus ^ and Poisson's equation for scalar poten- 
tial ([2]). The first of them describes the transverse field 
produced by the radiation incident on plasma, and the 
second one describes the longitudinal electrostatic field 
due to charge separation in plasma. 

The relativistic hydrodynamic equations for a plasma 
with uniform bulk electron temperature Tg and static 
ions can be written as 



dt 



div(7Vev) = 0. 



(4) 



l^ + (vV)p = -reViVe + - — +eVv?--[vxrotA] (5) 
at c at c 



7moV. 



(6) 



Here v, p are the mean velocity and pulse of the elec- 
tron component, respectively. This system may also be 
reduced in one-dimensional geometry. As a result, in the 
stationary case it reduces to two equations, the first of 
which is the law of conservation of a canonic electron 
pulse in transverse direction: 



c 



0. 



(7) 



The second equation is a simple balance of forces acting 
on the electrons: 



d7 



1 dn 
ndC 



0. 



(8) 



The first term here corresponds to the electrostatic force 
of charge separation in plasma, the second to the pon- 
deromotive force of the laser field, and the third term 
corresponds to the internal pressure forces arising due to 
temperature (/i = Te/rriQC^ is the parameter characteriz- 
ing electron temperature) . The continuity equation ^ is 
fulfilled automatically in this case, thus ensuring plasma 
quasineutrality. 

We first consider solutions to the system of equations 
(HI), ([2]) and ^ in infinite homogeneous plasma {ZiUi — 
1 throughout the space). In this case, it is convenient 
to integrate equation ([8]) and express the concentration 
through 



exp ■ 



'-7 + 1 



(9) 



where it is implied that we are interested only in localized 
solutions for which a— >0,<^^0,n^ lis fulfilled for 



C ±oo. In addition, s = also holds for localized 
solutions in infinite plasma. Substituting ^ into the 
two remaining equations yields the following fourth-order 
system of ordinary differential equations 



d^ 
d^ 



no 



: exp l a (10) 



no exp 1 



(11) 



Unfortunately, the above system has no analytical solu- 
tion, but we found one of its integrals that is written 
as 



1 /da\' 



1 



2 \dC J 



1 

2 VdC 



-fino exp ■ 



no4> - 



const. 



(12) 



This integral simplifies search for localized solutions to 
our problem. Indeed, let us set a problem of finding 
localized solutions symmetric relative to the C — plane. 
The relation 



da 







C=o 



is fulfilled for such solutions. And the relationship be- 
tween the values of vector and scalar potentials at the 
same point gives the integral (|12p : 



1 



a'{0)-nom- 



/ (/)(0)-v/l + a2(0) + l \ 
uno exp 1=0. 

V ^ / 

We now have one fitting parameter (for example, (/'(O)), 
by changing which and solving the system (|10p - pT|) we 
find a localized solution. Our studies showed that in such 
a formulation there may exist a family of solutions differ- 
ing by the number of field antinodes. Examples for the 
first and third modes are given in Fig[T] Note that, for 
inexplicable reason, we failed to construct modes higher 
than the first one for small values of parameter /i. This 
is most probably a result of numerical errors made in 
constructing the solution, but it may also have physical 
explanation. 

Besides, the integral allows us to show that there 
are no antisymmetric solutions within the framework of 
our system for which the relation 



= 



d(f) 
dC 



= 0. 



is fulfilled in the C = plane. At the same time, we 
constructed such solutions in a cold plasma (see, e.g., 
[l3|). This isn't surprising either physically or mathe- 
matically. Indeed, in terms of mathematics we increased 
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Figure 1: Examples of localized solutions of system (|10|l - 
(a) First mode, one field antinode. (b) Second mode, 
two field antinodes. (c) Third mode, three field antinodes. 



the order of the studied system, which may lead to appre- 
ciable changes in solution topology and phase plane upon 
the whole. In terms of physics, introduction of nonzero 
temperature results in appearance of electrons that can 
move from one field antinode to another during heat mo- 
tion, thus providing their interaction. This interaction 
most likely leads to the nonstationarity of antisymmetric 
distributions. 

Solutions constructed in an infinite plasma point to 
some principal features of the regime of interaction of a 
superintense laser pulse with overdense plasma. These 
are appreciable redistribution of electron concentration 
under the action of ponderomotive force and the related 
longitudinal electrostatic field approximately equal to the 
radiation electric field, as well as existence of structures 
with several electron layers separated by gaps almost de- 
void of electrons (fig. [T]d). 

Let us now pass over to the case of a bounded plasma 
layer that is more important for practical reasons. We 
will restrict our consideration to the one-sided irradia- 
tion of plasma. Then, it is convenient to construct so- 
lutions starting at some point outside plasma, from the 
unirradiated side rather far from the plasma so that the 
conditions n^(j)',a' << 1 (the primes mean differentia- 
tion along the coordinate) can be regarded to be met at 
this point. Note that a « -y/s, and the function 4>{C,) is 
generally unbounded (in the one-dimensional geometry 
the potential is known to behave logarithmically at large 
distances). This makes the system (flUl) - pT|) unhandy 
for constructing solutions. We remind the reader that 
this system was obtained from the basic equations ([2]), 
(l3|) and Q by omitting the variable n. By introducing 
a new variable e that denotes a longitudinal electrostatic 
field and is define by 



dC 

we obtain an equivalent system of equations. In this fash- 
ion we omit from the initial system of equations variable 
(j) for which boundary conditions are not clear and obtain 



the following system 
d^ _ 

de 



— no {Zini{C) - n) 



dVT 



dC 



dC 



(13) 
(14) 

(15) 



Here ni(C) is a given function specifying initial plasma 
distribution in space. We will further restrict our con- 
sideration to the case of a plane homogeneous layer of 
thickness L so that this function will have the form 







for C < 0, C > i^L/c, 
for < C < ujL/c. 



(16) 



.(0 - 



By analogy with ([TO]) - (HH) the system ^ - ([TS]) has 
an integral analogous to (I12p which, in the case under 
consideration, has the following form 



1 



2 (a) 



1 



—n^Zi I eni{C,)A(l ~ const. (17) 



Unfortunately, for rii(C) ^ the integral p7|) is not an 
algebraic expression and contains integration in explicit 
form. This restricts its applicability. However, it may be 
used in the region fully devoid of ions, where the last term 
in the left-hand side of the expression (fT7|) is identically 
equal to zero. 

As was mentioned above, we will construct solutions 
starting at some point C,st that is rather far from the 
plasma layer from its unirradiated side. Suppose that 
a = ^/s + a at this point and the inequalities OL(C,st) << 1, 
n{(st) « 1, e{(st) « 1 arc fulfilled. In addition, 
«i(Cst) = as this point is outside the plasma layer. 
Then, in the neighborhood of the point Cst we can lin- 
earize equations - ([T5|) in a standard manner to ob- 
tain 



dC2 
de 

dC 
dn 



-4a -|- riQn 



1 + s 



— = -non 



dC 



= — n e + 



da 



1 + s dC 



(18) 
(19) 
(20) 



Note that it would be wrong to neglect in the third equa- 
tion quantities of the form ne and nda/dC that are small 
compared to dn/dC because the resulting system would 
have only a trivial solution; hence we preserve these 
terms. As the system (fT8|) - (|20|) is still nonlinear, its 
solution is rather complicated, but assuming that the in- 
equality 



e >> 



da 



1 + sdC 



(21) 
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is met, equations (jT9l) and pO]) are separated from ((T8)) . 
which allows finding their solutions. Indeed, if we divide 
the second of these equations by n, differentiate it and 
substitute de/d^ from the first equation, we will obtain 
one second-order equation for n: 



1 /dn 



[dC 



non 



0. 



(22) 



Let us seek solution to this equation in the form 
and obtain 

^iCp {(3 - 1) (C - Co)''-' - f^cpHc - Co)"-'- 

-noC^iC-Cof =0. 

This equality must be identically fulfilled for all values of 
C; therefore, the conditions 



/3 = -2, 



no 



must be met. Finally the solution of equation ([22]) takes 
on the form 



n{0 



From equation (|19p we find 

2n 



e(C) 



C-Co' 



And from (fTSl) follows 



«(C) 



2(C - Co) V 1 + 



(23) 



(24) 



(25) 



We can now readily verify that the inequality ([2T|) that 
was supposed to be met does hold sufficiently far from 
the point C = Co- Note that all the integration constants 
in the derivation of (|23p - (P5|) were chosen so that the 
quantities a, e, n vanished at the infinity. 

Generally speaking, in the construction of the solution 
Co is an unknown quantity that determines a whole fam- 
ily of solutions. There arises a question which of these 
solutions is correct. The answer to this question is quite 
trivial. One must choose the solution at which the con- 
ditions of plasma quasineutrality upon the whole are ful- 
filled: 



+ 00 



(Z,n,(C)-n(C))dC = 0. 



(26) 



Actually, the problem of constructing a solution reduces 
to sorting all solutions by parameter Co and finding the 
one that would satisfy the condition (P^ . Moreover, 



there is no need to introduce Co in this case. Indeed, 
from - ([^5]) we can find the following relations 



(27) 
(28) 
(29) 



da 
dC 



2/ino 



1 



Thus, we have excluded parameter Co- 

Now the solution procedure will be the following. We 
choose a certain point Cst > i^L/c at which we set the 
value of variable e. The values of a, da/dC, and n at this 
point are found from the relations (P7)) - (|29p taking into 
account that a — ^/s-\-a. Thus, we have boundary condi- 
tions for the system (|13p - at the point Cst - We solve 
the system of equations and check whether the equality 
is fulfilled. If it is fulfilled, then it is the required 
solution; if the inequality is not fulfilled, we change the 
value of the variable e at the point Cst and repeat the pro- 
cedure. We proceed in this fashion until we find the right 
solution. In view of complexity of the system its solution 
was done numerically by manually sorting the value of 
e(Cst)- (There are some ideas of how this sorting might 
be automated, but they have not been implemented for 
lack of the barest necessity.) 

Now a few words about parameter s. As we have men- 
tioned above, it stands for the e.m. energy flux density 
and is equal to the intensity of radiation passing through 
the layer. In our problem it determines the amplitude of 
the wave incident on the layer. Thus, if we have to con- 
struct a function of the incident wave amplitude, then 
we must sort parameter s for each of its values using the 
procedure described above. 
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Figure 2: Example of solution of system ((13} - ([15} in a plasma 
layer of finite thickness. Parameters of the problem: = 0.01, 
no = 2, uL/c — 2, s — 0.1. The corresponding amplitude of 
incident radiation at = 0.85. 
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An example of solution constructed by the described 
algorithm is presented in fig. O The parameters of the 
problem were taken to be no = 2, s — 0.1, ujL/c = 2. 
As soon as the solution has been constructed one can 
find the amplitude of the incident pulse which, in a real 
situation, is a true parameter of the problem specified in 
the course of experiment. For this the values of a and 
da/dC must be known at an arbitrary point in vacuum 
from the irradiated side of the plasma layer. Indeed, at 
this point there are incident and reflected waves, hence 
we can write the following equality 



a(C)e 



(30) 



where ^{C) was introduced in equation ([T]) and and a^. 
denote complex amplitudes of the incident and reflected 
waves, respectively. By differentiating equation (|30p and 
omitting from the resulting system the quantity we 
obtain 



da . fdd ^ 



(31) 



Taking the modulus of this expression and bearing in 
mind that parameter s was introduced in the form (see 
0) 

odi^ 
dC 

we finally obtain an expression for the real amplitude of 
the wave incident on the layer: 



-I 1/2 



(32) 



It is also interesting to compare solutions with differ- 
ent values of electron temperature ji. If parameter s is 
fixed, then the value of the incident wave amplitude will 
be different for different values of fi. Therefore, such a 
comparison will not be physically correct. For a more 
correct comparison we fitted parameter s to each value 
of /i so as to ensure equal incident wave amplitudes in all 
the cases. The result is demonstrated in fig. El Note that 
the solutions obtained for /z = by the method described 
in [3] are also shown in this figure. The characteristic 
feature of this solution is enhancement of electron layer 
border blurring with increasing temperature, which is a 
natural and expected effect. Upon the whole, it should 
be noted that solutions of the system (fT^ - with 
< /i < 0.01 differ from the solutions with /i = ob- 
tained in the work ^14] only slightly and are comparable 
in terms of complexity of their construction. This en- 
ables us to conclude that both the approaches to find- 
ing steady-state distributions are equivalent. The only 
drawback of the solutions with nonzero temperature is 
divergence of the value of scalar potential at the infinity. 
However, this divergence is logarithmic and may be elim- 
inated manually by truncating solutions at some distance 
from the plasma layer. 
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Figure 3: Comparison of solutions of system (|13|l - (|15|l for 
different values of ^ but equal incident wave amplitudes. Dis- 
tributions of electron concentration are given on top, and of 
laser field at the bottom. The solution for /i = was con- 
structed by the method described in [14] . 



In [1J| it was shown that the so-called resonator-like 
structures in which electron layers separated merely by 
an ion gap play the role of walls may exist in a layer of 
finite thickness. Generally speaking, structures with an 
arbitrary number of electron layers separated by ion gaps 
may exist in a plasma layer. However, a special algorithm 
is needed for constructing solutions with N electron lay- 
ers as such solutions are discontinuous, with the number 
of solution discontinuities being proportional to N . This 
circumstance strongly complicates construction of such 
solutions, whereas under the action of temperature the 
solutions become smooth and all the family may be ob- 
tained identically by simply varying the values of s and 
e(Cst)- An example of solution with three electron layers 
one of which is outside the plasma layer is given in fig. U) 
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Figure 4: Example of solution of system (jlSp - (|15|l with three 
electron layers. Layer parameters: ^ = 0.01, no = 2, luL/c — 
3.4; incident wave amplitude ai = 2.1. 
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To conclude, we proposed a method for constructing 
steady-state structures in a plane plasma layer irradiated 
by superintense laser pulse. Introduction of uniform bulk 
electron temperature enables us to obtain smooth solu- 
tions. In spite of the fact that the problem is a model 
one, the obtained solutions demonstrate basic features of 
overdense plasma dynamics in the field of superintense 
waves, such as electron density redistribution under the 
action of ponderomotive forces, the presence of apprecia- 



ble electrostatic field in the longitudinal direction, and 
possible existence of resonator-like solutions with several 
electron layers separated by relatively rarefied gaps. The 
presented solutions may also be used for analysis of pos- 
sible ways of producing strong longitudinal fields in solid- 
state targets, as was discussed in, e.g., [l^[lB|- Besides, 
these solutions may be a handy tool for analysis of the 
phenomenon of relativistically induced transparency [l^ . 
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